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Abstract 

We compute the M theory corrections to the confining linear potential between 
a quark and an anti-quark in J\f = 1 Super Yang-Mills theory. We find a constant 
term, and a term exponentially small with characteristic length of Aq^, d . 

The potential in the MQCD setup that corresponds to softly broken M = 1 SYM 
is found to have a similar behavior. 
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1 Introduction 



During the last few years new light was shed on gauge theories via world volume field 
theories on D-branes stretched between two NS5 or D-branes, in IIA or IIB string theories 
HI (see also || and references therein). In particular, the brane configurations of D = 3+1, 
M = 1 Super Yang-Mills (SYM) and Super QCD were analyzed. The M theory versions 
of those theories, living on the M5 brane {|, |5|, ||, give a powerful new tool to deal with 
those theories (It should be stressed, though, that the M theory versions are not identical 
to the original theories, although they lie in the same universality class). In particular, 
the identification of the QCD string as an M2 membrane ending on the M5 |J, enables 
the calculation of the quark anti-quark potential via the MQCD string tension 0. The 
leading order of the potential was found to be linear in the separation / between the quark 
anti-quark pair, E oc I. This result which is valid, of course, only for large separations, 
(for I ^> Aq CD where the meson is "stretched"), is an indication of confinement. 

In this note we "measure" the potential with higher precision than that of 0. We 
derive the Nambu-Goto action that corresponds to an SU(N) gauge theories. In the case 
of SU (N = 2) we explicitly determine the sub-leading terms of the potential dependence 
on I, in the limit Aqcd < 1 (in M theory units). We find that the corrections to the 
linear term include a constant term, and a term exponentially small with characteristic 
length of Aqcd- We must remark, though, that because the MQCD string is not a BPS 
state, its tension is not protected when we vary the parameters of the theory. 

Similar nature of the sub-leading corrections to the potential are also detected in 
calculations involving MQCD setups that correspond to softly broken M = 1 SYM theory. 

Our work involves calculations similar to ones performed recently in the context of 



the AdS/CFT correspondence f8fl-fT8[. It turns out that the next to leading corrections 



found in our analysis are of the same nature of those found in a supergravity setup that 
corresponds to the pure four dimensional YM theory [13, Tl]. We also find no trace of a 
oc 1/7 Liischer [IJ3] term correction which is believed to exist in lattice simulations |2I| of 
YM theory. We argue that incorporating quantum fluctuations of the string world-sheet 
may induce such a term as was also the case in ref. [TJI . 



2 QQ potential in N = 1 MQCD 

The potential energy of a quark anti-quark pair translates in the MQCD language to 
the area of the M-theory membrane that corresponds to the QCD string. Therefore, 
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the extraction of the potential energy involves the construction of the MQCD setup 0, 
introducing a convenient parameterization of the membrane, and computing its area. We 
begin by reviewing the MQCD setup of 0. 



2.1 The setup 

The MQCD description of M = 1 SQCD is given in terms of a five-brane (M5) embedded 
in R 10 x S 1 = R x i? 4 x Y , where the 11-th compact dimension S 1 is a circle of radius R. 
The M5 has the form i? 4 xE, where £ is a Riemann surface in the space Y which is a three 
dimensional complex space with coordinates v = £4 + ix$ , w = xj + ixs , t = e~ s l R where 
s = xq + iRxio ■ The setup does not extend in the additional non compact dimension, 
and we can express this as the requirement xg = 0. 

In the limit R — > the M5 reduces to the well known type IIA brane configuration 
of figure (HD |ZT| . In this configuration the flavor degrees of freedom are associated with 



the semi-infinite D4 branes stretched to the right of the NS5 brane. An alternative setup 
includes semi-infinite D4 branes stretched both to the right of the NS5 brane and to the 
left of the NS5' brane. This later configuration is considered in section 2.4. Although we 
have quarks in our setup they are very heavy and we take for E the Super Yang-Mills 
SU(N) curve |, §, § 



v N 



v-w = ( (1) 

C N 'H (2) 



where ( is a complex constant. We measure distances in M-theory units, namely, in units 
where the Planck length is set to be l p — 1. We see that for each t there are N distinct 
(v,w) points, differing by powers of the phase e 2m l N ( and therefore there are iV branches 
of the M5 brane. 

The MQCD M5 is characterized by two dimensionful parameters, R and (, whereas 
the M = 1 SYM has only a single scale Aqcd- Matching the string and domain wall 
tensions of MQCD with the corresponding ones of the large N SYM requires the following 
assignments ||: 

ICI ~ N 2 A* QCD (4) 

The dependence of R and ( upon Aqcd is dictated from dimensional analysis. As we shall 
deal with the small N regime, and are not interested in numerical factors in the Aqcd 
language, the N dependence of (§),(|D wm n °t be used by us. 
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Quarks and W-bosons are described in MQCD by membranes wrapped around the 11- 
th dimension. This leads to open strings of type IIA in the limit R — > 0. The boundaries 
of such membranes must be curves on E with zero winding- number around the 11-th 



dimension |P2" |. Such membranes may have the topology of a disk or of a cylinder. It 
231 that membranes with the cylinder topology correspond to vector 



was shown 

supermultyplets, while those with the disk topology correspond to chiral ones. In |7], B3| 
it was demonstrated that in our setup we cannot describe a membrane representing a 



are confined. 
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Figure 1: The M = 1 branes configuration 
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Figure 2: The QCD string for N=6. this is the v plane for a constant xq 

The construction we are interested in is of a "meson" , that is, of a quark and an anti- 
quark connected by a QCD string. The membrane corresponding to the QCD string is 
stretching in one direction along R A , while in the other direction it interpolates between 
two adjacent branches of the M5 brane (see figure (0)). This membrane will not vary 
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in the x\o direction, so we might take x±o = 0, and have s = x§. The length of the 
interpolating segment for a certain s, which is the width of the membrane, is (as we will 
see shortly) 

2 v^C sin(7r/iV) y/ cosh(2s / RN) (5) 

(From here on ( stands for \(\). The membrane corresponding to a fundamental type 
IIA string ("quark") transforms smoothly into a QCD string and back into a fundamental 
string. The membrane's area is the total energy of the "meson" in M-theory units. This 
area tends to infinity as the quarks' masses approach infinity, so we need to subtract 
those masses. We consider only the part of the membrane that behaves as a 'QCD string 
membrane' rather than 'quark membrane' (fundamental string in the type IIA string 
theory). We first characterize the point of the transition - this is the point where the 
tension of the QCD string is the same as the type IIA string tension: 

2 v /2C v / cosh(2s/i?A0 = 2ttR (6) 

Once we make all our calculations up to this value of s, the quarks, whose masses we 
wish to subtract, are simply QCD strings with constant coordinates in R A and s, going 
from the maximal value we found (which we will eventually take to infinity) to zero. 

2.2 The area of the membrane 

The QCD string is associated with an M2 stretched in the Y space, as well as in one 
spatial dimension of R 4 , that we take as x = 23. For each (x, s), the membrane's slice is 
a straight line segment in (v, w) with both ends on E. We have the freedom to rotate the 
v plane so that this segment will vary only in x$, and likewise for w and x$. Therefore, 
we can write a parameterization of the membrane as 

P(x,p) = (x 3 , x 4 , x 5 , x 6 , x 7 , x 8 ) 

= (x , ^/(e- s{x)/RN cos(7r/iV) , p^fC,e~ s{x)/RN sm^/N) , (7) 
s(x), ^/(e s{x)/RN co S (ir/N), -py^e s ^ RN sin(7r/iV)) 

with — 1 < p < 1. For p = ±1, P is on E and the distance between these two points gives 
the string's width (||). 

The function s(x) is determined by minimizing the area of the membrane. The width 

(in the (v , w) space) of the QCD string decreases as s approaches 0, so the membrane has 

2 Actually, in this parameterization, = —Q N / 2 e~ s / R ^ so it corresponds to x\q — n. This, however, 
has no physical significance, and we shall continue to treat s and xq interchangeably. 
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an "incentive" to achieve small values of s. On the other hand, of course, this involves a 
large detour in s. The best compromise is achieved, as we shall see, for a minimal value 
of s, being so = s(0), exponentially small in I. 

Taking the membrane's area as an action, we can write a Lagrangian which is an area 
element in the form 



£ 



— I A I — 

dp J \ dx 



det[h a f- 



where h a p = d a P ■ dpP, and a and (3 correspond to p and x. Inserting the explicit 
expression of (|7|) one finds 



£ = V^Csin(7r/iV)v'6 + ap 2 



(9) 



where 



2( 



R 2 N'- 



sin 2 \n/N)s 



,2 



COS 



h(2s/RN) (1 + s' 2 + cos 2 (n/N) cosh(2s / RN)s' 2 



(10) 



'111 



integrating over p we get : 



£ = ^/2Csin(7r/iV) (Va + 6+ ^= (log(Va + 6 + v 7 ^) - log(V a + b - v 7 ^))^ 



12) 



We look at the limit Aqcd *C 1 (in M-theory units 0), so by (QfJ), for not too big an 



c 



R 2 N 2 



< 1 



(13) 



and therefore: 



£ = 2 v/2~C sin(vr /N) ^/cosh(2s /i?iV) y 1 + s' 2 + cos 2 (vr/A^) cosh(2s/#iV)s /2 



(14) 



which has a very simple form for 577(2): 



£ = 2 v /2C v / cosh(s/7?) v 7 ! + s' 2 



(15) 



3 in those units, the eleven-dimensional Planck length is set to one, l p = 1. The dimcnsionful condition 
is Aqcd < X - 
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This form has a simple interpretation. It corresponds to the area ignoring the change 
in the (v, w) coordinates of P(x, p = ±1), as a function of x, and retaining only its (x, s) 
change. Vl + s' 2 is then the length of the line element of the QCD string in the (x, s) 
plane, and 2y / 2£y / cosh(s/i?) is the width of the membrane in the transverse (v, w) plane. 
This width, then, provides naturally the metric, which is scalar diagonal. 



However, although the change in x^ and xj is suppressed by the factor (0), we cannot 
ignore it if it tends to infinity with s. The only case in which we can ignore this change 
is when 0, that is, the case N = 2. Hence, this interpretation is valid only for 

N = 2, where, indeed, X4 = 27 = 0. 

The Hamiltonian of the N = 2 action is 

H _ -2y^C V /cOsh( g /i?) 



vi + s' 



We shall first treat the case in which the quark and the anti-quark are on the same 
NS5 brane (in the II A picture), that is, they both reside in the s — > +00 limit. Due to 
the symmetry of the problem we have s'(0) = 0, and defining s = s (0) we can write (as 
the Lagrangian does not depend explicitly on x): 



H(s, s') = 2 ^CVcosh(s/R) = _ 2 ^- Wcosh{so/R) = H{ Q) 
Vl + s' 2 

This differential equation can be written as: 



which gives : 



s'(x) = ±v^cosh(s/-R)/ cosh(s /-R) — 1 



x = ±l - (19) 
y cosh (s/R)/ cosh ( sq / R) — 1 



In particular, since by (H,|l3D we have s(l/2) ~ log(R/ ^ 1> we can write: 



1 1 = r y/ C Osh(s /R) ds 

2 J so A/cosh(s/-R) — cosh(so/-R) 
Inserting the solution we found into the action, we see that the energy of the QCD 
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string is: 



2 



C dx 

«('/2) ^ x 

C—ds 

ds 



SO 

2a/2C 



cosh (s/i?) ds 



s A/cosh(s/i?) — cosh(s /-R) 

1 /"°° 

= 2 v /2C v /cosh(s / J R)(-0 + 2 V / 2C / V cos H s /R) ~ cosh(s /R) ds 

2 Js 

from which we must subtract the (infinite) mass of the quark : 

POO 

M q = 2^/2( y/cosh(s/R) ds 
Jo 



(21) 



(22) 



changing the integration variable to y — cosh(s/i?) — cosh(s /-R) and z = cosh(s/R), 
and defining e = cosh(s /-R) — 1, we finally get: 

dy 



I = 2Ry/(l + e) 



o y/y{y + e){y + 2 + e) 



/l+e 



dZy/z 

Vz^T 



where 



J(e) = f 
Jo 



dy 



y/V - Vv + 1 + 1 

^{y + e){y + 2 + e) 



(23) 
(24) 

(25) 



1/2 



Figure 3: A characteristic shape of s(x). The asymptotes represent "bare" quarks 



2.3 The potential dependence on the QQ separation 

The distance I as a function of e is shown in the appendix to take the following form 



1= ^-v / 21og£ + 5v / 21og2- -^j=e\oge + 0(e)^ R 



(26) 



The explicit dependence of the energy E on e, which is also computed in the appendix, 
reads 

E = 2^2? I (l + \< - \* + (e >) + 2 W 0(-« + ^ , log i - ^ S log ± + G(e ^ 
with k = -J(0). 

Inserting the value of £ from the inverse of (p6|), we finally find the following expression 
for E as a function of / 

E = 2^/2(1 (l - 2{R/1)k + 9Qe~ 2l/V2 ~ R + 0((R/l)e~ 2l/V2R )^j (28) 

It is easy to see how the corrections to the linear potential depend on Aq C d, the 
parameter of M = 1 SYM theory, by translating the MQCD parameters according to (0,^). 
As a function of I, one correction is constant, while the other is essentially exponential 
in —l/R ~ —Aqcd I- In fact, given the dependence on I, the dependence on Aqcd is 
determined, since Aqcd is the only energy scale of the original theory. 

Note that the potential does not include a Luscher term, namely, a term of the form 
—c/l, where c is a positive universal constant. Such a term is missing also in the Wilson 
line derived from the Ads/CFT correspondence H]. Since the traces of such a term were 
observed in lattice simulations (and in phenomenological fits to the heavy mesons 
spectra), the authors of ref. |14[ suggested that this may indicate that there is a phase 



transition. The proposed scenario is that the Luscher term exists in the weak coupling 
phase but is missing in the strong coupling phase. It seems to us that the origin of the 
discrepancy is different and, in fact, is hidden in the work that led to the discovery of the 
c/l term [TJ. The 

source of this term is the contribution of the Gaussian fluctuations 
around the classical solution. In the present work, as well as in those that follow the 
Ads/CFT duality, the computation was based solely on the classical solution and the 
quantum fluctuations have not been incorporated. It is plausible that the evaluation 
of the fluctuations around the classical configuration in the stringy calculations of both 
MQCD and the AdS/CFT duality will produce a term of the form c/l. 



2.4 QQ potential of the symmetric MQCD setup 

We proceed now to discuss the case in which the quark and the anti-quark are, asymp- 
totically, at s = ±oo. In the IIA picture, they are on the two NS5 branes, as opposed to 
the setup discussed so far, in which they are on the same NS5 brane. 

The previous case gave rise to an even function s(x), whereas now we get an odd 
function. We expect the energy for this case to be higher than the previous one. The 
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reason is that if we reflect, for negative x, the graph of the odd s(x), we get an even 
membrane with the same area. The solution found in the previous sections for the even 
case has, of course, a smaller area. (The reflected membrane would not be in equilibrium 
for x = 0, s(x) = 0). 

The Lagrangian and Hamiltonian are equal to those of the previous case, but, as 
s(0) = 0, the solution is now specified by the value of r = s'(0), and 

= -2VT ( ^m = _ 2M i = 

VI + s' V 1 + r o 

If we define for this case e = r 2 and substitute y = cosh(s/i?) — 1, we get 

I = 2R [°° dy (30) 

and the energy, after subtracting the quark mass, is 



E = 2 



i +m ,i) r iy ^HE^^vF±i) (31) 



Using methods similar to those used in the even case, we get for the energy a stronger 
dependence on e. This time we have 

E = 2^2(1 (l-\e + 0(5 2 ) + 2(R/1)(-k - ^= e log ± + ^=e + 0(e 2 loge))^ ^ 

We can, therefore, write / only to the first order in e. This happens to be identical to 
the even case, so we have for e 

e = 32e- l / VlR + 0(e~ 2l /^ R ) (33) 

Inserting this into the expression of the energy, we have for this case: 

E = 2^2C, I (\ - 2(R/1)k + 2(R/l)^=e- l ^ R + 0(e- 2l ^ R )^j (34) 

This potential has the same type of behavior as the one with the asymmetric MQCD 
setup. The corrections to the string tension term are a constant and an exponential term. 
However, it is clear that for small enough e (That is, large enough /) this result is larger 
than the one we got for the even s(x) case. 
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3 QQ potential in softly broken Af = 1 MQCD 



We proceed now to analyze the QQ potential in a MQCD setup that corresponds to softly 
broken M = 1 SYM. In the field theory picture this is achieved by giving mass to the 
gaugino and in the context of brane configurations this is done by rotating one of the NS 
branes in a non-holomorphic way (e.g. rotation in the (2:4, £7) plane). The notations that 
we are using are those of || where the brane is parameterized by two four dimensional 
complex null vectors p and q, 

p 2 =q 2 =0 (35) 

in the following form 

(x 4 ,x 5l x 7 ,x 8 ) = Re(p\ + q\^ 1 ) (36) 
x 6 = -RNcRehgX (37) 
x w = -iVImlogA (38) 

where A is a complex number and c is a real number given by 

P 2 AT 2 

-P"<f+^-(l-c 2 ) = (39) 
In this notation we choose the M = 1 curve to be 

p= VC(l,-i,0,0) «f= Vc(0,0,l,-z) (40) 

Using the rotation suggested in |2R EBj we keep p the same as for Af = 1 and choose 



q 



\/C(sina, i sin a, cos a, — i cos a) (41) 



For simplicity, we shall concentrate now on £77(2) YM, namely, considering the brane 
at a constant s = e~ X6 (that is, A = re l6 = s 1 / 2 ^ e je ) we find the cross-section with the v 
plane to be a circle of radius y/C {r + sin a / r), and the cross section with the w plane a 
circle of radius ^/Ccosa / r. The distance between two points of equal Xio is therefore 



2 v /2C v / cosh(s/i?c) + sin a = 2^/2( cosh(s /Rc) + p ■ q (42) 

Using the arguments given above for the supersymmetric £77(2) theory we can write 
the Lagrangian of this theory, defining /x = sin a, as 

C = 2 v /2Cv / cosh(s/i?c) + fiVl + s' 2 (43) 
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and the expressions for the length and the energy of the QCD string 

dy 



PC 

I = 2R^{l + n + e) / 

Jo 



E = 2^/2(^1 + fx + el + A^2(R J{e) - 4y/2(R j 



v /y(y + e)(y + 2 + e) 



(44) 
(45) 



where 



J{e)=l dy ^-^ + l+ » + j (46) 
lo v /(y + e)(y + 2 + e) 



and we used (according to ([TBI)) 



c 



^l-p-q^R 2 = y/1 - ^(/R 2 ^ 1 (47) 



Using the same methods as in the supersymmetric case we find 

log(£/32) = -7kvr^( 1 - 32( ^4»- ,/(OT, ) +0(e " /(OT,) 
£/32 = e - mm + i| < _L= ( _i_ _ m 

+ o(e~ 2l/{V2RVTJ ^ ) ) 
and for E(l(e),e, /x) 

^ /— r , / / £ £ 2 -R, e , e e 2 , e .\ 

£ = h/l + W+ ; ; ^ + 2^(-Ku + = log = log— ) 

V ^ V 2v^+^ 8(1 + /i) 3 / 2 Z l 2^ 32 32^ 8 32^ 

+ 0(Re 2 ) + 0(fe 3 ) (49) 

where — k m is the zeroth order term of J. We finally have for the energy in the softly 
broken SYM theory 



E = 2 ^2C I (VT+Ji - 2{R/1)k (U + /2 (-5 - y. + 8^/1+^)6-^^^) 

+ 0(( J R//)e- 2Z/(v ^ /? ^ I) ) N ) (50) 



Using the relation of the rotation parameter a to the gaugino mass given in 26 1 we 



get in our notations the plausible result rri\ = tana, and the pure YM theory is given by 
a = 7r/2 and \i — 1 . In this case the energy takes the form 



E = 2^/2~C I (V2 - 2(R/1)k 1 + 16(8 - 3V2)e~ l/R + 0{{R/l)e- l,R )\ (51) 

Comparing this result to the one of the M = 1 case we see that it has the same princi- 
pal functional dependence on I. The string tension is modified, however, and so is the 
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exponential correction, since the screening length is R rather than R/ v2 as in the super- 
symmetric case. As the equations (Hf|) give the relation between R, ( and Aqcd only upto 
dimensionless constants, the significance of those modifications in field theory is unclear. 

An interesting special case of this calculation is when jj, — — 1, when the curve spans 
only the Xq and v coordinates and is singular at |A| = 1 (or Xq = 0), where the corre- 
sponding QCD string's width is zero. The behavior we find for this case is very similar to 
the one described in Jy] (A/" = 4 at finite temperature). Once we regulate the integrals 
(by taking the upper limit as s max instead of oo) we find that for / ^> y/ s max R the quarks 
are free (the favorite configuration is of two "bare" quarks with the QCD string located 
at x 6 =0) while for I <C V s maxR we still have linear confinement as the leading order 
behavior. 



4 Summary and Discussion 

Quark anti-quark potentials have been extracted recently from computations of Wilson 
loops in various supergravity models which are conjectured duals of supersymmetric gauge 



theories in the large N limit |j8[-||18|. It is thus interesting to compare the calculations in 
the MQCD framework to those of the gauge/gravity dual models. In both cases, in spite of 
the fact that the starting points seem very different, the calculations reduce eventually to 
computing the area of a bosonic string world-sheet via a Nambu-Goto action. Moreover, 
using a gauge a = x where a is a world-sheet coordinate and x is a coordinate along which 
the quark anti-quark pair is situated, the MQCD and gravity calculations are special cases 



of an integral of the following form |28 



S = / dx^f%s(x))+g%s(x))(d xS y (52) 



where s is a coordinate external to the R 1,3 space-time, s = Xq in the MQCD case and 
the fifth coordinate of the AdS§ space (or of its relatives). In fact, in the approach of 
Polyakov |2"Tfl , s is associated with the Liouville coordinate. It is straightforward to realize 
that f 2 (s), g 2 (s) are both cosh(s/R) in SU(2) MQCD, (s/R Ad s) 4 , 1 in the Ads/CFT case 
§, {s/R Ads )\l - (s T /s) 4 ), 1 in finite temperature [|IJ and (s/R Ads ) 4 , (1 - (st/s) 4 )" 1 for 



the "pure YM" 4 dimensional model |13fl . Naturally we would like to compare the quark 
potentials extracted from the stringy computations. In particular the MQCD result for 
broken supersymmetry ( |5TD can be compared with the precise computation of pure YM 



theory in four dimensions (eqn. (22) of |T4|)- 111 both cases the corrections to the linear 



potential do not include a c/l Liischer term, but only an exponential term of the form 
le~ bl , where b is proportional to Aqcd- We believe that a c/l term will show up once 
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quantum fluctuations are incorporated as was the case in |Tj| Unlike the MQCD result, 



there is no constant term in the expression derived from the supergravity setup [I4||, but 
this may follow from a different subtraction. 

In this work on precision "measurements" of the quark potential from MQCD, we 
found two types of corrections to the linear potential between external (very heavy) quark 
and anti-quark in SU(2), M — 1, SYM. The first is a constant correction proportional to 
R (the radius of the 11-th dimension) and the second is essentially exponential in —AqcdI- 
We considered two configurations - one with both quarks on the same NS5 brane, and 
the other with the quarks on the two opposite NS5 branes. Although these configurations 
are supposed to describe the same field theory, we found the corrections are different. We 
looked at theories where supersymmetry is softly broken and found that the behavior is 
qualitatively the same - a linear potential with a constant and exponential corrections. 
For more general SU(N) (S)YM one can write analogous equations to ( |23| , |24| ). However, 
as is clear from (|23|), for N ^ 2 the expressions will be much more cumbersome. 
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5 Appendix - The explicit computation of the asymp- 
totics of / and E 

5.1 Asymptotics of / 

The computation involves extracting the asymptotic behavior of 1(e), and substituting its 
inverse in the asymptotic expression of E(s). The leading behavior of 1(e) reads 

1= ^-v / 21oge + 5v / 21og2-^=eloge + 0(£)^ J R (53) 

In order to prove that, we first show that 

t(\ P dy 1 51og2 1 

his) = / = = = log£H = — =eloge + 0(e) 

Jo ^y(y + e)(2 + y + e) V2 8^ (g4) 

This behavior can be calculated in a few steps. First we look at an integral with 
the same divergence for e ~ as h(e), for which we can find a primitive function and 
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therefore evaluate its value as a series in e : 

ji\ r ^ d V !, , 3\og2 1 

/2(£) "X vxv + m + 2)' = "vi l0££ + - wf £ l0££ + Q(£) (55) 

All that remains is to evaluate the contribution of the difference l\ — Ii. As both 
integrals have the same divergence for e pa the difference is continuous for e = and its 
leading contribution is given by setting e = : 

A/(0) = h(0) - J 2 (0) = r^^±^f? = V21og2 (56) 

io y( 2 + y) 

The e log £ term of the difference is found by looking at the derivative with respect to 



d - AI(e) = [°° dy~ £ ^ + 2y + + e + - (2 + y)(2 + 2y - y/2y/2 + e + y) 



e : 



de w J » 2Vy(2 + y)(2 + e + y)s/*(e + y)«/a 

Again, we shall look for an integral with the same divergence for e ~ 0, and having a 
primitive function. This time we have: 

l°° dy ,— ~stn) ^ = ^log£ + 0(l) (58) 

Jo y ^y(e + y) 3 / 2 (2 + y) 3 8v/2 1 ' 

(As before, checking that the difference of the integrals does not diverge when we set 
e = provides a proof that we found the correct diverging term). Integrating the result 
we got with respect to e, we find that the difference's contribution is £^oge, and show 
that Ii{s) is as in (]54). Inserting this expression in (|23|) 

/ = 2 J R V /(TT7) /°° ^ (59) 

Jo ^y(y + e)(y + 2 + e) 

we finally find (0). All that is left is to invert this expression, namely 

log(e/32) = --^ (l - 12e-'/^) + 0(e~ l ^ R ) (60) 

e/32 = e- l l^ 2R + ^{l/R)e- 2l ^ m + 0{e- 2l ^ m ) (61) 
V2 

5.2 Asymptotics of £" 

The expression for E(l(e),e) was given in (p4[) 

£ = 2 v^Cv/T+i I + 4v^Ci? J(e) - 4v/2Ci? / 1+£ -^£= (62) 

Ji yz 2 - 1 
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It is straightforward to convert the first term into 

2 v /2Cv / TT7/= ^l + ~£ + 0( £ 2 )J 2^/2(1 (63) 

For the third term we use z ~ 1 within the whole integration range so that we can 
express the integrand as a series in z — 1: 

= v^-v^+^ 3/2 + 0(r>/ 2 ) (64) 

The second term of (|23|) has a finite negative contribution for e = 0. We define Q 

k = - J(0) » 1.63029 (65) 

and so this finite contribution is 

i^/2( R J(0) = -2k ■ 2y^2( R (66) 

To find the sub-leading order behavior in e we write the integral J(e) as a sum of 
three integrals 

J(e) ee J 3 ( £ ) + J t ( e ) - J 2 ( £ ) 

rfy v^-v^ 7 ^ + Z 1 dy ^ ( 67 ) 

i V( y + e )(2 + y+ i) 7o V(l/ + e)(2 + y + e) 



/o ^(y + e)(2 + y + £) 
Js(e) is an analytic function at e = 0, so it is possible to write it as a series in e 



/oo 
dy 



Vviy + 2) 

dy ( - + (1 + ^7 CT ' ) e + 0(s*) (68) 

^2VTTV^2+^) Z/(2 + l/)v / R2T^) ' 



re can also be given a closed form in terms of elliptic functions. 
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The behavior of J\{s) and ^(s) can be found in a similar way to the one we used to 
find I: we find a function with the same divergence, subtract, differentiate with respect 
to e and repeat the process . Using this method we find: 

( 1 log (2 + VQ) 3 , \ 1 

Jl(£) = Cl+£ (7! + yH-7! log2 j + v! £log£ 

' e 2 \oge + 0(e 2 ) (69) 



Me) = C 2 - V2 ■ V~e + - -^e 3/2 + 0(e 2 ) (70) 

Summing up all the contributions to the energy we have indeed shown Q27| ) 
B = 2V2C( (l + i £ - I? + Otf + 2( B /,)(-K + ^ e log ^ - ^ .og^ + 0(^)) 
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